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INTROUUCT  ION 


This  report  presents  .in  evaluation  of  the  proposed  kort  nozzle  for  tlu  I'.i 
WYTM  Cutter.  ItV'onsiilers  the  shroud  structure,  the  propeller  hydrodynamics, 
t  tie  ice  1  outline,  on  l  lie  propeller  blade,  a  possible  i  ee  deflect  ion  strut  ,  and 
overall  design  of  the  kort  nozzle.  The  shroud  structure  is  considered  first  in 
Section  l,  wucre  it  is  found  to  be  sufficiently  strong.  -Sect-ton  2  goes  into 
the  propeller  hydrodynamics,  while  Section  3  considers'  t he  remaining  topics. 

One  serious  problem  witli  the  present  design  is  that  there  is  no  access  to 
the  propeller.  Unless  the  port,  starboard,  or  both  sides  of  the  nozzle  is  remo¬ 
vable,  the  propeller  would  have  to  be  installed  and  the  shroud  built  around  it. 
The  screw  would  be  inaccessable  for  repairs.  Suggestions  for  how  this  problem 
could  be  solved  are  put  forward, in  Section  1. 

The  thrust  and  efficiency  of  the  propeller  together  with  the  shroud  meet 
requirements.  However,  the  local  velocity  distributions  are  not  hydrodynarai- 
cally  the  best  and  result  in  cavitation  on  portions  of  the  back  and  face  of  the 
blades.  Slight  modi f ica t ion s  to  the  sections  should  alleviate  this  situation;  a 
very  efficient  subcavitating  propeller  should  result.  Details  may  be  found  in 
Section  2  and  Appendix  C. 

Unless  the  propeller  blades  are  to  be  steel,  they  are  not  strong  enough  to 
withstand  extreme  ice  Loading.  Details  of  the  ice  loading  are  given, in  Section 
3  and  Appendix  D.  Recommendations  4  and  5  in  Section  4  address  this  problem. 

Supports  for  the  nozzle  are  also  discussed  in  Se-’inn  T.  .  '.nose  should  be 
faired  in  to  prevent  adverse  hydrodynamic  drag,  which  reduces  the  total  net 
thrust. 

Recommendations  may  be  found  in  Section  4. 
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1.  Till'.  SHROUD  STRUCTURE 


To  predict  the  ice  loading  to  which  tin1  nozzle  will  he  s  ub  jec  ted  is  not 
possible;  however,  the  loading  will  not  exceed  the  crush  in;;  strength  of  the  ie,. 
Gr.il  l  (1977)  suggests  using  the  crushing  strength  of  the  ice  times  the  pro¬ 
jected  area  as  the  upper  bound  of  the  compressive  force  that  can  be  t  ransr.ii  tted 
from  the  ice  to  a  structure.  He  also  suggests  the  use  of  an  overall  factor  of 
safety  ranging  from  1.5  to  2.5.  Garivick  and  Lloyd  (1970)  suggest  the  use  of 
3ik  p.s.i.  for  the  crushing  strength  of  the  ice.  Since  there  is  no  way  to  anti 
cipate  the  total  area  against  which  the  ice  will  act  nor  its  location,  we  will 
use  Graft's  criterion  for  the  loading.  The  nozzle  will  be  subjected  to  300 
p.s.i.  acting  horizontiallv  as  shown  in  Figure  1-1,  since  this  direction  would 
be  the  most  critical  Loading.  It  is  assumed  that  the  nozzle  will  be  rigidly 
supported  at  the  hull;  this  assumption  will  be  discussed  in  detail  in  Section  3 


Each  half  of  the  nozzle,  the  port  side  and  starboard  side,  will  be  analyzed 
as  an  arch  with  fixed  ends.  Figure  1-1  shows  both  the  applied  forces  and  the 
reaction  forces  that  act  on  each  half.  Roark  (1954)  gives  the  following 
expressions  for  the  reactions: 

a)  the  horizontal  reaciton  H 


h)  the  vertical  reaction  V 


V  =  WR  (I  -  I  +  2±)  (1.2) 

_ 6 _ 4 _ 8 _ 

2_  (H_  -  1 )"  -  31  +  2  -  all 
!i  2  4  4 

c)  the  moment 

M  =  WjT  -  2HR  (1.3) 

4  n~ 

where  R  js  the  radius  of  the  neutral  axis,  w  the  loading  per  linear  foot, 
d  =  1/ AR  f  and  A  is  the  cross  sectional  area. 

Section  dimensions  and  other  properties  of  the  shroud  are  tabulated  in 
Appendix  A.  Th e  chord  c  of  the  nozzle  is  4.282  ft;  hence,  the  loading  per 
linear  foot  is: 

w  =  300  p.s.i.  *  144  in2  ft  *  4.282  ft 
w  =  185,000  lbs/ft 

The  neutral  axis  of  the  arch  is  the  circle  through  the  centroid  of  the  cross 
section;  its  r_adius  R  is  measured  from  the  center  of  the  propeller  shaft  and  has 
the  value  R  =  y  =  4.6^  ft.  The  section  area  A  =  2.194  ft  ,  and  its  moment  of 
inertia  I  =  0.0961  ft  . 

Substituting  these  values  into  eauations  (1.1)  thru  (1.3)  yields 

a  =  0.961  =  0.0020 

7.194"*  '4767^ 

H  =  185,000  lbs/ft  *  4.67  ft 
=  864,000  lbs 

V  =  864,000  lbs  *  (-0.0825)/ (-0.1503) 

=  474,000  lbs 

M  =  (864,000/4  -  2  *  474,000/n )  *  4.67  ft  lb. 

=  -  401,000  ft  lb 

The  maximum  compressive  stress  obc  due  to  bending  is  given  by: 

°bc  =  il  max  "  R) 

I 

=  401,000  ft,  lbs  {  5.34  -  4.67)  ft 


°bc 


=  2,796,000  p.s.f. 
=  19,400  p.s  .i . 
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There  is  an  additional  compressive  stress  oac  due  to  the  axial  load  H  at  the 
support . 

°ac  =  H/A 

=  864,000  lbs  /  2.19  ft2 
=  304,000  n.s.f. 

°ac  =  2,740  P-5-1 • 

The  maximum  compressive  stress  is: 

°c  =  °ac  =  Jhc 

=  2740  p.s.i.  +  19400  p.s.i. 

=  22,100  p.s.i  . 

At  the  fixed  support  the  maximum  tensile  stress  due  to  bendina  is 

aj,t  =  M  max  (y]  -  R) 

I 

=  401,000  (4.67  -  4.37)  p.s.f. 

0961” 

=  1,251,000  p.s.f. 
ofot  =  8690  p.s.i. 

Since  there  is  an  axial  compressive  stress,  the  maximum  tensile  stress  is 
qiven  by: 

°t  =  °bt  '  °ac 

=  8,690  p.s.i.  -  2,740  p.s.i 
at  =  5,950  p.s.i. 

The  shear  stress  t  is  qiven  by: 

T  =  v/A  =  474,000  lbs  /  2.194  ft2 
=  216,000  p.s.f. 

=  1,500  p.s.i . 

The  internal  bendinq  moment  around  the  arch  may  be  expressed  as: 

M(  e)  =  -  M  -  VR  sine  +  HR  (1-cose)  -  wR2  (1-COS0)2 

~T~ 

=  -  M  -  VR  sine  +  HR  (l-cos9)  -  HR  (1-cose)2 
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Where  M  posit  ive  indicates  inner  side  (  lower  surface)  in  tension. 


(  9)  =  -  HR  cosQ  +  VR  sinO  -  VR(l-cosu)  sinu 

=  -  HP  cosO  +  VP  cose  sinu 

The  roots  of  M'  ( o)  =  n  are  0  =  n  and 

sine  =  HP  =  H  =  474,000  =  0.5486 
VR  V  864,000 


or 


0  =  0,  33.27°,  n  -  33.27° 

The  moments  at  these  points  have  the  values 
M( 0)  =  206,000  ft  lb 

M  (33.27°)  =  M  (rr  -  33.27°)  =  -206,000  ft  lb 

The  level  of  stress  for  the  arch  is  small,  so  that  a  low  carbon  steel  can  be 
used  for  its  construction.  This  has  the  additional  advantaqe  that  low  carbon 
steels  are  a  touqh  material  with  good  resistance  to  impact  loadino.  The 
nozzle  will  be  subjected  to  impact  loadings  from  ice  at  low  temperature. 

With  the  present  nozzle  design  there  is  no  access  to  the  propel ler.  Ti  nee 
the  stress  level  within  the  shroud  is  low,  the  port  or  both  +u-  port  and  star¬ 
board  sides  could  be  made  removable  to  provide  the  neces^ar .y  access.  One 
approach  would  be  to  make  the  nozzle  a  built  up  structure  with  a  cross  section 
1  ike  Figure  1.2. 


The  structure  would  consist  of  a  leading  edge  beam  sculptured  to  the  profile  of 
the  section  and  a  trail ina  edqe  beam  also  sculptured. 

llnper  and  lower  surfaces  between  the  beams  would  be  made  from  plates. 

There  would  be  ribs  at  approximately  45°  intervals.  The  outer  and  inner  surface 
plates  would  be  made  removable.  The  leading  and  trailing  edge  beams  would  be 
made  removable  near  the  ship  hull  and  the  lower  support.  One  approach  to  these 
connections  would  be  to  have  oermenent  short,  beams  attached  to  the  hull  and  the 
lower  support;  the  leadina  edqe  and  trail inq  edqe  beams  would  bolt  to  these 
fixed  structures.  If  this  proves  impossible,  the  spacing  between  the  two  beams 
mioht  he  made  1 arge  pnouqh  for  the  prop  to  fit  throuqh.  The  removable  nozzle  is 
however  preferable. 
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2.  PPOPFLLFP  HYDRODYNAMICS 


Fuller  (1981b)  provided  the  desiqn  criterion  for  the  screw.  At  a  ship  speed 
Vs  =  6  knots  the  propeller  is  to  turn  at  N  =  270  rpn  at  a  shaft  horsepower  SHP  = 
2400  hp.  TTie  wake  reduction  factor  is  to  be  w  =  0.177;  this  relates  to  a  velo¬ 
city  at  r/R  =  0.6.  The  axial  velocity  into  the  propeller  is  then  qiveri  as: 

Va  =  (1-w)  Vs  =  8.340  ft/sec 

Details  of  the  hydrodynamic  force  calculations  may  be  found  in  Appendix  C. 

At  the  above  conditions  the  thrust  developed  by  the  screw  is  aiveri  in  Table  C-8 
as: 


T  =  33,100  lbs 

and  the  corresponding  toraue  is: 

0  =  46,800  ft  lbs. 

These  values  are  for  a  water  density  p=1.94  lbs.  sec2/ft.  At  this  toraue 
the  delivered  horsepower  would  be: 

DHP  =  2  it  nO  =  2,354  Hp 
550 

Where  n  =  4.5  r.p.s.  is  the  rotational  shaft  speed. 

s  =  DHP  =  0.981 
SHP 


Which  is  right,  where  it  should  be. 

The  thrust  coefficient  for  the  screw  alone  is: 

Kt  =  _  T.  =  0.162, 

Tscrew  p7f^- 

and  the  torque  coefficient  is: 

Kq  =  _0.=  0.026 
pn2  D5 

Van  Manen  and  Superiana  (1959)  determine  the  thrust  constant  of  the  nozle 

KT„02z1e  '  °-110' 

consequently,  the  total  thrust  constant  is: 

^screw  +  ^nozzle  "  °-272 

The  total  thrust  is: 

T  =  pn2  D4  KT  =  59,300  lbs, 

and  the  correspondi no  thrust  horsepower  THP  is: 

IMP  =  T_Va  =  899  Hp 
550' 


Overall  efficiency  for  the  screw  no77le  combination  is: 

=  T  HP  =  0.38? 

'°  D  HP 

At  258  rpm  for  a  delivered  horsepower,  DHP  =  ?350  Hp,  the  bollard  Dull  due 
to  the  proneller  is: 

TSCrew  =  37,400  lbs 

The  thrust  coefficient  for  the  screw  is: 


Addina  this  to  Kt  ,  vields  a  total  thrust  coefficient 
1  no77l  e  • 

Kj  =  0.309 
for  a  total  thrust 

T  =  58,129  lbs 

This  compares  with  the  present  bollard  pull  of  56,000  lbs  at  245  rpm.  The 
shroud  has  reducted  the  bollard  pull  of  the  screw  by  a  third  without  lowerinq 
the  total  bol 1 ard  pul  1. 

The  velocity  distribution  is  obtained  through  use  of  the  Goldstein  (IS 52) 
approximation  to  Theodoresen's  mapping  technique.  Velocity  dicfr  luutions  for 
r/R  =  0.5,  0.7,  and  0.9  are  shown  in  Figures  2-1  through  2-3.  Details  of  this 
method  can  be  found  in  Appendix  C.  Although  the  velocity  distribution  for  the 
section  at  r/R  =  0.5  at  a  Cl  =  1,  figure  2-4,  is  quite  normal;  in  contrast,  at 
the  design  value  Cl  =  0.146  the  lower  surface  velocities  exceeds  those  on  the 
upper  surface  for  the  first  15%  of  the  chord.  At  r/R  =  0.7  this  condition 
holds  over  the  first  2.5%  and  at  r/R  =  0.9  over  the  first  10%  but  not  as  drasti¬ 
cally  as  at  r/R  =  0.5.  One  other  bad  feature  of  the  velocity  distribution  is 
its  behavior  over  the  last  25%  of  the  upper  surface  in  the  vicinity  of  r/R  = 

0.7.  This  erratic  behavior  is  still  evident  for  the  same  section  at  Cl  =  1. 

The  propeller  cavitates  under  simmer  operating  conditions.  Details  of  the 
cavitation  calucul  ations  are  in  Appendix  C.  There  is  cavitation  on  the  back  of 
the  blade  at  r/R  =  0.7.  There  is  cavitation  along  the  lower  leading  edqe  from 
an  approximate  r/R  =  0.5  to  the  tip.  Both  of  these  conditions  are  due  to  poor 
velocity  distribution  rather  than  to  an  excessive  blade  loading.  It  should  be 
possible  to  operate  the  propeller  cavity  free,  after  minor  modifications  to  the 
section  shape. 
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FIGURE  2-3  VELOCITY  DISTRIBUTION  r/R=0.9 


).  OTHER  STRUG  i'l-RAl.  AND  HYDRODYNAMIC  CONSIDERATIONS 


Stresses  due  to  ice  crushing  by  the  blade  are  tabulated  i  a  Table  IW>  and 
were  derived  in  Appendix  D.  These  stresses  exceed  the  design  working  stresses 
for  both  mangancese  bronze  and  nickel  nlmninun  bronze  screws  recommended  bv 
O’ lirien  (I9h2).  In  fact  the  stress  levels  would  actually  call  for  a  steel 
screw.  Some  might  argue  that  the  crushing  stress  of  ice  is  too  severe  a  loading 
condition.  consider  the  results  for  r/R  =  0.8;  the  stress  8800  p.s.i.  repre¬ 
sents  a  chunck  of  ice  being  forced  against  the  outboard  1 o. 2  inches  of  tin- 
blade.  True,  the  block  would  have  to  be  s’  wide,  but  this  is  plausible.  The 
bl.ule  could  not  support  an  even  smaller  block  pushing  normal  to  its  face-;  for 
at  r/K  -  0.8  the  appropriate  1  is  7  in  ,  the  maximum  distance  is  0.039  feet  or 
0.468  inches,  and  using  an  a llow  stress  of  6500  p.s.i.,  the  maximum  moment  would 
be 


^a  l  lowed 


6300  p . s . i.  *  7  in  4 
0.468  in 


97.222  in  lb 


A  lo  inch  thick  by  10  inch  wide  piece  of  ice,  acting  with  a  3  inch  moment  arm, 
would  mearly  have  to  exert  194  p.s.i.  to  achieve  this  level  of  stress. 

A  strut  stretched  between  tiie  hull  and  the  shroud  as  shown  in  Figure  3-1 
could  be  used  to  deflect  ice  away  from  the  blades.  It  would  not  however  protect 
the  blades  from  ice  coming  in  astern  during  backing  operations.  An  XACA  0018 
section  would  be  sufficiently  strong  to  withstand  ice  loading,  as  t he  predicted 
maximum  stress  would  be  7780  p.s.i.,  see  Appendix  i).  It  is  not  possible  to  pre¬ 
dict  the  direction  of  flow  into  the  strut;  a  visual  study  in  a  water  channel 
would  be  advisable  to  check  out  the  hydrodynamics  of  the  -trut  and  i ^  interac¬ 
tion  with  the  propeller. 


FIGURE  3-1  ICE  DEFLECTION  STRUT 
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For  hydrodynami c  reasons  the  upper  and  lower  supports  for  the  nozzle  should 
he  faired  into  the  ship  hul  1.  Strucutrally  they  should  be  made  an  inteqral  part 
of  the  hull,  so  that  t  h*\y  art.  as  rigid  supports  for  the  nozzle.  Figure  3-7 
shows  a  '  ut  through  t  tie  upper  support  (shaded  portion)  to  show  the  fairing.  The 
width  of  tiie  support  should  he  determined  by  the  location  of  its  lower  outside 
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1  SUPPORT 
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Figure  3-2  ATHWART  SHIP  CUT  THROUGH  UPPER  SUPPORT 


edge;  this  edge  should  be  located  so  that  it  just  provides  clearance  enough  tr 
remove  the  propeller.  This  is  approximately  the  2'  6"  buttocks  plane.  Figure 
3-3  shows  a  typical  fore  and  aft  section.  Two  major  features  :rc  the  slots  into 
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Figure  3-3  FORE  AND  AFT  SECTION 


which  the  shroud  leading  edge  beam  and  its  trailing  edge  beam  fit.  These  beams 
should  bolt  to  the  support.  It  may  be  useful  to  have  an  access  plate  at  mid 
chord  to  access  these  bolts;  this  view  shows  this  plate.  Note  also  that  the 
support  continues  to  the  rudder  support  housing  and  should  be  faired  into  this 
structure . 
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Figure  3— ♦  shows  a  stern  view  of  a  suction  through  the  lower  support . 


ourstb£ 

SUKPAUF 


Figure  3-4  LOWER  SUPPORT  SECTION  STERN  VIEW 
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4.  RECOMMENDATIONS 


1.  The  nozzle  structure  should  be  made  removable  to  provide  access  to  the 
propel ler. 

?.  Either  the  pitch  angle,  the  slope  of  the  mean  line,  or  both  should  be 
adjusted  so  the  sections  operate  closer  to  optional  CL  and  the  unusually  high 
velocities  on  the  face  near  the  leading  edge  are  reduced. 

3.  Section  thickness  distributions  should  be  investigated  to  determine  if 
they  cause  the  cavitation  on  the  back  of  the  blade.  Velocity  distributions  are 
notproperly  distributed  and  the  section  shape  needs  modification  to  alleviate 
the  situation. 

4.  If  the  present  cutter  propellers  have  not  received  any  significant  ice 
damage,  these  propellers  should  be  analyzed  to  determine  the  maximum  stresses 
they  will  support  at  different  r/R  values  for  a  safety  factor  of  two.  These 
values  could  then  be  used  as  ice  load  design  stresses  for  the  propeller  under 
consideration  in  this  report. 

5.  Either  strengthen  the  blades  to  support  the  loading  described  in 
Appendix  D  or  design  to  the  stresses  described  in  recommendation  4. 
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APPENDIX  A  NOZZLE  ANALYSIS 

The  shroud  dimensions  in  this  report  are  from  the  140  WYTM  Cutter 
Propel,  ler  Data  which  are  stored  on  the  Naval  Academy  Time  Sharing  System 
(NATS).  This  data  was  digitized  from  the  United  States  Coast  Ouard  Drawing 
(Fuller  1981a).  The  section  data  presented  in  Table  A-l  is  good  to  0.005  ft 
i.e.  SO. 06  inches.  Both  yu  and  y^  are  scaled  data;  whereas,  the  fore  and  aft 
coordinate  X.  is  a  direct  reading  in  inches  from  the  drawing.  x  is  the  scaled 
version  of  X,  where  the  scale  factor  1  ' / 1 . 5 ”  is  used. 


TABLE  A-l  SHROUD  SECTION 


Upper  Surface 


Lower  Surface 


X(in) 

x(f  t) 

Yu(ft) 

X(in) 

x  (ft)  j 

YL(ft) 

3. 176 

2.117 

5.1  14 

3.176 

1 

2.117 

5.114 

3.125 

2.083 

5.  180 

3.038 

2.025 

4.906 

3.025 

2.017 

5.233 

2.621 

1.747 

4.650 

2.925 

1.530 

5.240 

2.183 

1.455 

4.481 

2.236 

1.491 

5.150 

1.584 

1.056 

4.371 

-0.022 

-0.0  15 

4.918 

0.774 

0.516 

4.314 

-0.022 

-0.015 

4.925 

-0.106 

-0.070 

4.316 

-0.092 

-0.061  S 

4.898 

-1.536 

-1.024 

4.330 

-1.790 

-1.193 

4.730 

-2.586 

-1.724 

4.372 

-3. 1 18 

-2.125 

4.590 

-3.106 

-2.701 

4.402 

-3.228 

-2.1 52 

4.549 

-3.207 

-2.138 

4.436 

-3.247 

-2.165 

4.482 

-3.247 

-2.165 

4.482 

TABLE  A-2  presents  this  section  data  as  a  function  of  the  nondimensional 
distance 

x1  =  (x|_f  -  x)/c 

where  xLy  is  the  x-coordinate  of  the  leading  edge.  From  the  table  xLp  =  2.117 
ft;  c  is  the  section  chord  length  and  is  given  as 

c  =  XLE  -  XTE 

T.E.  denotes  trailing  edge;  from  the  table  xyp  =  -2.165  ft.  Hence, 
c  =  2.117  -  (-2.165)  =  4.282  ft. 

In  the  table  n.t.  denotes  nose  to  tail  line,  i.e.  the  line  between  the  leading 
edge  and  the  trailing  edqe.  Yc  is  the  camber  line  approximation: 

Yc  =  0.5  (Yu  +  Y-| )  -  n.t. 

The  upper  surface  coordinate  Yu  and  lower  surface  coordinate  Y1  are  linear 
interpulations  of  the  data  in  Table  A-l.  The  maximum  thickness  of  the  shroud  is 
0.712  ft  from  the  tabulated  values  of  Yu-Yj ,  and  this  mimimum  thickness  occurs 
at  25^  of  the  chord.  Maximum  camber  also  occurs  at  25%  of  the  chord  and  has  the 
absolute  value  0.230  ft. 


G  cm  r\  C  m  vc 
C  h  ro  a'  co  kc 
O  O  O  G 


r—4  C  fs"»  C\i 

G  h  ro  w  fv 

C\j  C\J  C\J  C\J  C\J 


od  ro  r — » 

h  cc  isir. 

C\-l  r—4  r—4  f — I  *— 4 


>, 
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^  C  C\-  O";  CO  IT)  O'  f— I  VD  CO  H  ai  fO  IT)  in  H 

H  Oi  C  ^  CC  HNCvitJlN  ^  ro  H  CC 

f— iOOOO^OC’  CO  N  N  l£)  VO  VO  VC  VO  uo  VO 


xn  vx>  tx>  u~> 


^  cv  o  o 
vo  in  c  ro 

LX,  LT,  IX)  IX. 


O'  r— i  f-H  C\J 

»-i  a  o  cc 

IX)  IX) 


Section  properties  for  the  shroud  are  calculated  by  dividing  the  shroud  sec¬ 
tion  into  elements  ax  =0.05  units  wide,  as  shown  in  Figure  A-l. 


Figure  A-l  0 IVI SION  OF  SHROUD 

Since  yu-y|=0  at  both  the  leading  and  trailing  edges,  the  section’s  total  area  A 
in  terms  of  the  trapazoidal  rule  is: 

A  =  0.05c  E  (yu-yi ) 

By  definition  and  the  trapazoidal  rule  the  location  ~y  of  the  section  centroid 
from  the  axis  of  the  propeller  shaft  is: 

y  =  0.05c  l  yi  (yu-yi ) i/A » 

where  yi=(yu+yi ) /2.  Since  the  centroid  of  each  element  is  not  located  at  the 
same  distance  _from  the  propeller  axis  as  the  centroid  of  the  section  but  is 
additional  yi-y  away,  it  is  useful  to  use  the  parallel  axis  theorem.  When  this 
is  done,  the  moment  of  inertia  T  around  the  centroidal  axis  for  the  shroud  sec- 
t ion  i s  given  by: 

T-  o.05c  (I(yu-yi)(yi-y)2+ni 

where  Ti  =  (yu-yi)3/12. 
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TABLE  A- 3  SECTION  PROPFRTIFS 
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APPENDIX  R:  PROPELLFR  BLADE  GEOMETRY 


The  propeller  dimensions  presented  in  this  report  are  obtained  from  the  140 
WYTM  Cutter  K  3-75  Propeller  Data  stored  on  NATS.  This  data  was  diqitinized 
from  United  States  Coast  Guard  Drawina  (Fuller  1981a).  The  section  data  pre¬ 
sented  in  Tables  B1  throuqh  B7  is  good  to  _+  0.005  ft.  which  is  the  same  as 
_+  0.06  in.  The  propeller  radius  P  =  4.25  ft. 

TARLE  B1  niGITIZED  SECTION  DATA  r/R  =  0.4 


Upper 

Surface 

Lower 

Surface 

x  (ft) 

vtl  (ft) 

x  (ft) 

Y]  (ft) 

1.851 

0.198 

1.851 

0.198 

1.850 

0.21* 

1.835 

0.180 

1.835 

0.226 

1.769 

0.182 

1.768 

0.256 

1.597 

0.176 

1.656 

0.304 

1.426 

0.177 

1.481 

0.359 

1.297 

0.174 

1.323 

0.404 

1.108 

0.166 

1 .113 

0.450 

0.859 

0.166 

0.851 

0.499 

0.622 

0.170 

0.633 

0.519 

0.386 

0.168 

0.393 

0.529 

0.152 

0.170 

0.187 

0.532 

-0.005 

0.173 

0.001 

0.534 

-0.269 

0.165 

-0.00? 

0.532 

-0.566 

0.175 

-0.049 

0.527 

-0.860 

0.170 

-0.?38 

0.516 

-0.991 

0.171 

-0.421 

0.501 

-1.132 

0.174 

-0.744 

0.464 

-1.334 

0.183 

-0.965 

0.418 

-1.514 

0.187 

-1.783 

0.345 

-1.689 

0.194 

-1.482 

0.290 

-1.750 

0.196 

-1.603 

0.259 

-1.753 

0.200 

-1.726 

0.220 

-1.756 

0.204 

-1.754 

0.210 

-1.756 

0.204 
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TABLE  B-2  DIGITIZED  SECTION  DATA  r/R  =  0.5 


Upper  Surface 


x  (ft) 


Lower  Surface 


Yu  (ft) 

x  (ft) 

0.193 

1.988 

0.199 

1.984 

0.209 

1.966 

0.223 

1.932 

0.258 

1.831 

0.295 

1.686 

0.332 

1.448 

0.381 

1.248 

0.424 

1.047 

0.449 

0.686 

0.461 

0.162 

0.464 

0.001 

0.468 

-1.181 

0.467 

-0.408 

0.457 

-0.638 

0.443 

-0.911 

0.430 

-1.207 

0.397 

-1.396 

0.369 

-1.618 

0.319 

-1.810 

0.260 

-1.935 

0.2  23 

0.206 

0.196 

-1.941 

Yi  (ft) 


93 

75 


0.167 
0.165 
0.162 
0.161 
0.165 
0.165 
0.166 
0.165 
0.169 
0.167 
0.170 
0.166 
0.172 
0.176 
0.174 
0.172 
0.172 
0  191 
0.182 
0.182 


TABLE  8-3  DIGITIZED  SECTION  DATA  r/R  =  0.6 
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Upper 

Surface 

Lower 

Surface 

X  (ft) 

Yu  (ft) 

x  (ft) 

Y]  (ft) 

1.990 

0.191 

1.990 

0.191 

1.980 

0.205 

1.978 

0.175 

1.946 

0.217 

1.938 

0.165 

1.823 

0.246 

1.852 

1.701 

0.266 

1.722 

0.161  ! 

1.556 

0.289 

1.485 

0.159 

1.217 

0.328 

1.265 

0.161 

0.688 

0.378 

0.956 

0.167 

0.372 

0.396 

0.612 

0.166 

0.001 

0.40  j 

0.222 

0.165 

0.001 

0.399 

-0.001 

0.167 

-0.134 

0.402 

-0.241 

0.167 

-0.416 

0.392 

-0.620 

0.172 

-0.616 

0.386 

-0.983 

0.172 

-0.807 

0.370 

-1.322 

0.173 

-1.037 

0.354 

-1.630 

0.172 

-1.315 

0.322 

-1.892 

0.170 

-1.490 

0.300 

-2.052 

0.171 

-1.712 

0.266 

-2.098 

0.173 

-1.855 

0.245 

-2.108 

0.180 

-2.002 

0.214 

-2.072 

0.201 

-2.100 

0.193 

-2.108 

0.180 

TARLE  B-4  DIGITIZED  SECTION  DATA  r/R  =  0.7 


Upper  Surface 


Lower  Surface 


X  (ft) 

Yu  (ft) 

x  (ft) 

Y 1  (ft) 

2.242 

0.181 

2.242 

0.181 

2.238 

0.191 

2.240 

0.169 

2.196 

0.195 

2.228 

0.163 

1.982 

0.222 

2.150 

0.162 

1.766 

0.242 

2.023 

0.163 

1.494 

0.274 

1.872 

0.166 

0.971 

0.314 

1.678 

0.163 

0.467 

0.330 

1.404 

0.164 

0.001 

0.333 

1.158  i 

0.166 

-0.001 

0.334 

0.640 

0.168 

-0.404 

0.331 

0.309 

0.170 

-0.680 

0.325 

0.116 

0.168 

-0.952 

0.315 

-0.092 

0.173 

-1.174 

0.301 

-0.382 

0.170 

-1.478 

0.276 

-0.707 

0.178 

-1.659 

0.261 

-0.989 

0.175 

-1.804 

0.248 

-1.326 

0.175 

-2.052 

0.215 

-1.570 

0.171 

-2.167 

0.202 

-1.850 

0.173 

-2.227 

0.192 

-2.046 

0.171 

-2.233 

0.184 

-2.179 

0.174 

-2.227 

0.178 

i _ 

-2.233 

0.184 
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E 


fnmi\ 


Upper  Surface 


Lower  Surface 


x  (ft) 

P.321 
2.315 
2.277 
2.162 
2.018 
1.920 
1.706 
1.399 
0.937 
0.508 
0.002 
-0.002 
-0.069 
-0.282 
-0.603 
-1.002 
-1.395 
-1.568 
-1.756 
-1  .925 
-2.256 
-2.314 
-2.324 


vu  (ft) 

0.1  78 
0.188 
0.190 
0.201 
0.208 
0.215 
0.227 
0.243 
0.266 
0.282 
0.284 
0.286 
0.282 
0.282 
0.279 
0.272 
0.2  53 
0.241 
0.232 
0.220 
0.198 
0.181 
0.181 


x  (ft) 

2.321 

2.309 

2.243 

2.168 

2.037 

1.880 

1.680 

1.452 

1.192 

0.745 

0.4R0 

-0.003 

-0.320 

-0.691 

-1.090 

-1.368 

-1.576 

-1.931 

-2.122 

-2.282 

-2.322 

-3.324 


Y!  (ft) 

0.1  78 
0.164 
0.159 
0.163 
0.162 
0.164 
0.165 
0.165 
0.169 
0.169 
0.169 
0.171 
0.167 
0.178 
0.179 
0.172 
0.173 
0.170 
0.1  72 
0.175 
0.175 
0.181 


TABLE  B-7  DIGITIZED  SECTION  DATA  r/R  =  1.0 


Upper 

•  Surface 

Lower 

Surface 

x  (ft) 

Yu  (ft) 

x  (ft)  j 

Y]  (ft) 

2.388 

0.172 

2.388 

0.172 

2.382 

0.178 

2.386 

0.160 

2.353 

0.182 

2.376 

0.160 

2.228 

0.191 

2.337 

0.166 

2.061 

0.196 

2.259 

0.161 

1.918 

0.199 

2.166 

0.167 

1.526 

0.214 

2.035 

0.164 

1.202 

0.217 

1.829 

0.164 

0.682 

0.223 

1.287 

0.164 

-0.002 

0.224 

0.786 

0.168 

-0.002 

0.227 

0.208 

0.169 

-0.129 

0.229 

0.001 

0.169 

-0.373 

0.223 

-0.511 

0.171 

-0.588 

0.227 

-1.074 

0.177 

-0.826 

0.227 

-1.524 

0.173 

-1.177 

0.220 

-2.007 

0.177 

-1.586 

0.213 

-2.177 

0.174 

-1.826 

0.203 

-2.300 

0.171 

-1.979 

0.202 

-2.376 

0.170 

-2.163 

0.194 

-2.380 

n.  173 

-2.300 

0.193 

-2 .374 

0.186 

-2.380 

0.178 

The  expanded  blade  data,  presented  in  Table  8-8,  was  derived  from  Tables 
R-l  through  R-7.  The  leading  edne  L,  E.  is  the  largest  positive  value  of  X,  and 
the  trail  ina  edge  T,  E.  is  the  largest  negative  value.  Camber  c  is  the  distance 
between  the  leading  and  trailing  edges.  Maximum  thickness  for  this  propeller 
blade  is  at  X  =  0;  the  value  in  the  table  is  computed  using  the  maximum  value  of 
Yu  and  the  value  of  Y]  nearest  X  =  0.  All  values  for  r/R  =  0.3  were  taken  from 
the  drawing  by  Fuller  (1981a). 


TABLE  P-8  EXPANDED  BLADE  GEOMETRY 


R  =  4. 

35  ft 

r/P 

•"(ft) 

L.E.  (ft) 

T.  E.  (ft) 

c  (ft) 

t  (ft) 

0.3 

1.275 

1.68 

-1.52 

3.20 

0.424 

0.4 

1.700 

1.851 

-1.756 

3.607 

0.361 

0.5 

2.125 

1.988 

-1.941 

3.929 

0.301 

0.6 

2.550 

1.990 

-2.108 

4.098 

0.236 

0.7 

2.975 

2.242 

-2.233 

4.475 

0.161 

0.8 

3.400 

2.321 

-2.324 

4.645 

0.115 

0.9 

3.825 

2.373 

-3.374 

4.474 

0.073 

1.0 

4.250 

2.388 

-2.380 

4.768 

0.055 

It  is  convenient  to  have  all  section  data  as  a  function  of  +'.e  nondimen- 
sional  distance 

xLr  -  x 

X  •  =  - 

C 

from  the  leading  edge  for  the  purpose  of  computing  the  pressure  on  each  winq 
section.  This  data  is  tabulated  in  Tables  B-10  thru  B-15.  In  the  tables  n.t. 
denotes  the  nose  to  tail  line  or  chord  line;  it  is  the  line  between  the  leading 
edge  and  the  trailing  edge.  yc  is  the  camber  line  approximation 


Yc  =  °-5  (yu  +  yi)  - 

The  upper  surface  coordinate  yu  and  the  lower  surface  coordinate  yj  are  1  inear 
i nterpul ations  of  the  respective  data  in  Tables  B-l  through  P-7. 
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TABLE  B-9  SECTION  DIMENSIONS  AS  FUNCTIONS  OF  NON  DIMENSIONAL  DISTANCF  FROM  L.E .  r/R 


CO  C  C  CO  C\J  •£! 

o  c:  r^.  vc 

oc;  c  co  vc  cc 

CM  CM  VC  CO 

CT;  r— <  C\J  PO  LT  UD 

O'i  H  (V  ^ 

^  Ul  IT.  ^  CO 

CO  CM  o  co 

«-H  CM  C\J  CM  C\J  CM 

cm  co  co  co  co 

CO  CO  CO  CO  co 

CO  CO  co  CM 

00  0*  <T»  CG  OG  VC  rHVCr- 1  OCh^-MO 

co  lt>  cji  cm  ld  o^-coom  c\j  co  co  cm  .-h 

f-H  CM  CM  CM  CO  CO  LT;  U~>  LO  LT,  If)  LT/ 


.215  0.361  0.178  0.183  0.270 
.395  0.314  0.184  0.130  0.249 
.576  0.266  0.190  0.076  0.228 
.756  0.204  0.204  0.000  0.204 
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TABLE  B-ll  SECTION  DIMENSIONS  AS  FUNCTIONS  OF  NON  DIMENSIONAL  DISTANCE  FROM  L.E.  r/R 


TABLE  B- 12  SECTION  DIMENSIONS  AS  FUNCTIONS  OF  NON  DIMENSIONAL  DISTANCE  FROM  L.E.  r/R  =  0.7 
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TABLE  R- 1 3  SECTION  DIMENSIONS  AS  FUNCTIONS  OF  NON  DIMFNSIONAL  DISTANCE  FROM  L.E.  r/R 


TABLE  B-14  SECTION  DIMENSIONS  AS  FUNCTIONS  OF  NON  DIMENSIONAL  DISTANCE  FROM  L.L.  r/P  =  0.9 
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TABLE  8-15  PiTCH  AND  PITCH  ANGLES 


Standard  K  3.75 


3.70 

0.4619 

4.40 

0.4119 

4.92  1 

0.3635 

5.51 

0.3439 

5.95 

0.3183 

6.38 

0.2986 

7.21 

0.30O 0 

8. 36 

0.3131 

24.79 

22.38 

20.23 

18.98 

17.66 

16.63 

16.70 

17.38 


2.2  6 
2.89 
3.53 
4.3? 
5.19 
8.2? 
7.25 
8.44 


JiO_WYTM  Cutter 


0.2821 

0.2706 

0.2644 

0.2696 

0.2777 

0.2912 

0.3017 

0.3161 


1 5.75 

15.14 

14.81 

15.09 

15.52 

16.23 

16.79 

17.54 


APPENDIX  C:  PROPELLER  FORCES 


The  approximate  solution  for  the  angle  of  zero  lift  given  by  the  Pankhurst 
formula  (abbott  1959) 

ao  2/C  'i  A  yc 

is  used  to  obtain  the  theoretical  zero  lift  angle.  The  values  of  A  are  tabu¬ 
lated  in  Table  C-l,  while  those  of  ag  are  tabulated  in  Table  C-2. 

TABLE  C-l  PANKHURST  CONSTANTS 


X1 

A 

0.000 

1.45 

0.025 

2.11 

0.050 

1.56 

0.1 

2.41 

0.2 

2.94 

0.3 

2.88 

0.4 

3.13 

0.5 

3.67 

0.6 

4.69 

0.7 

6.72 

0.8 

11.75 

0.9 

21.72 

0.95 

99.85 

1.00 

-104.90 

Burril  (1944)  presents  a  correction  factor  Ka  with  which  to  correct  the 
theoretical  zero  lift  angle  to  the  one  found  experimentally.  The  corrected 
angle  a0  actual  is  a^so  tabulated  in  Table  C-2.  Thickness  t  and  chord  are  from 
Table  B-8;  whereas,  the  percent  of  chord  at  which  maximum  camber  occurs  is  taken 
from  Tables  B-9  through  B- 15.  The  factor  Ka0  is  taken  from  Burris  (1944). 
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The  ship  speed  Vs  =  <>  knots,  the  KPM  =  N  =  2 7 (J  r|jn,  and  the  wake  reduction 

factor  W  =  U -  17  7  td  an  r/R  =  0.6,  and  the  shaft  horsepower  SHP  =  2500  lip  are  from 

United  States  Coast  Guard  notes  (Fuller  1981b).  Witli  the  usual  design  process 
each  blade  section  would  have  seen  an  axial  velocity  Va ,  where: 

V.,  =  (l  -  W)  vs 

=  0.823  *  6  kts 

=  4.9  38  kts 

V.j  =  8.3  40  ft/ sec 

However,  as  discussed  in  the  body  of  the  text  this  axial  velocity  is  multiplied 
by  a  factor  Va  X/Va  to  account  for  the  experimentally  determined  velocity  profile 
in  the  nozzle  at  the  location  of  the  screw  (  Van  Manen  1962).  From  Appendix  A 
the  nozzle  chord  length  =  4.282  ft.  and  the  screw  diameter  D  is  8.50  feet. 

The  values  of  Vax/Va  in  Table  C-3  are  for  a  1/D  ratio  of  0.50  and  were  obtained 

from  Van  Manen  1962.  The  rotational  velocity  Vrot  is  given  by 

Vrot  “  r  n 

where  n  =  270  rpm/60  sec/min  =  4.5  rps.  Hie  velocity  V  is  the  incident  velocity 
that  the  blade  section  sees; 

V  =  (Vax2  +  Vrot2)V2 

it  comes  at  the  blade  from  an  angle  i3  where 
8  =  arc  tan  (Vax/Vrot) 
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TARIF  C-3  INCIDENT  VELOCITIES  AND  RESULTANT  ANCLE  H 


r/R 

^ax^a(  ft/sec) 

vax(ft/ sec] 

- 1 

vrot  (ft/ sec) 

— 

V  (ft/sec) 

- f 

S(deg) 

0.3 

0.76 

6.34 

36.05 

36.60 

9.97 

0.4 

0.80 

6.67 

48.07 

48.53 

7.90 

0.5 

0.84 

7.01 

60.08 

60.49 

6.65 

0.6 

0.90 

7.51 

72.10 

72.49 

5.94 

0.7 

0.96 

8.01 

84.12 

84.50 

5.44 

0.8 

1.06 

8.84 

96.13 

96.54 

5.25 

0.9 

1.20 

10.01 

308.15 

108.61 

5.29 

1.0 

1.35 

11.26 

120.17 

120.69 

5.35 

Figure  C-l  shows  the  relationship  between  these  velocities,  3,  and  the 
pitch  angle  9. 


FIGURE  C-l  INCIDENT  VELOCITIES 


For  computational  purposes  it  is  useful  to  work  with  an  angle  $0  defined  by 


❖o  =  *  +  ant  +  aoact 

whpre 

apt  =  arc  tan  [(YLF  -  YTE)/C] 

The  angle  oq  is  the  pitch  anqle  of  the  nose  to  tail  line;  it  is  tabulated  in 
Table  C-4. 


TABLE  C-4  PITCH  ANGLE  <j>0 


r/R 

❖(deg) 

C(ft) 

YLE(ft) 

YtE(ft) 

ant(deg) 

a0( deg) 

*0 

0.4 

15.14 

3.607 

0.198 

0.204 

-0.10 

3.85 

18.8 

0.5 

14.B1 

3.929 

0.193 

0.192 

0.01 

3.33 

18.1 

0.6 

15.09 

4.098 

0.191 

0.180 

0.15 

2.63 

17.8 

0.7 

15.5? 

4.425 

0.181 

0.184 

-0.04 

1.55 

]  7.0 

0.8 

16.23 

4.645 

0.178 

0.18] 

-0.04 

1.03 

17.2 

0.9 

16.79 

4.747 

0.168 

0.180 

-0.14 

0.84 

17.4 

1.0 

17.54 

4.768 

0.172 

0.178 

-0.07 

0.65 

18.1 

The  pitch  angle  4>  is  from  Table  B-15;  the  chord  c  is  from  Table  B-8;  YLfr  and 
YyE  are  from  Tables  B-9  through  R- 1 5 ;  and  aoflct  is  from  Table  C-2. 
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Our  first  imperical  correction  is  to  account  for  the  cascade  effect.  This 
appears  as  a  correction  to  the  pitch  angle  <j>0;  though,  it  is  an  actual  correc¬ 
tion  to  the  anqle  of  zero  lift.  The  corrected  pitch  anale  <f>0  or  is  given  by 

$0  cascade  =  <j>0  -  Kg«0  *  a0  act  (C.l) 

where  the  correction  factor  Kga0  is  a  function  of  $j,  an  anqle  which  is  ini¬ 
tially  undetermined  and  o  .  The  values  of  Kqa0  are  found  in  Burril  1944.  If  V 
denotes  the  resultant  velocity  and  w  the  induced  downwash,  then  Figure  C-3  shows 
the  relationship  between  the  various  velocities  and  Bj. 

a  =  3  c  /  (it  D  r/R)  (C.2) 


Vro  t 


FIGURE  C-2  SECTION  INDUCED  VELOCITY 


The  determination  of  <|>0  cascade  is  part  of  the  total  interative  process  to 
determine  the  induced  downwash  w.  In  order  to  start  the  process  3]  is  initially 
assumed  to  be  equal  to  3.  4>0  cascade  is  calculated  based  on  this  value,  and  the 

induced  anqle  of  attack  aj  can  be  taken  as  4>0  cascade  -  3.  Figure  C.4  shows  the 
relationship  between  the  induced  angle  of  attack,  the  angle  of  zero  lift  a0,  and 
the  anqle  of  attack  of  the  blade  a. 


FIGURE  C .3  INDUCED  ANGLE  DF  ATTACK 

Rather  than  solve  for  the  downwash  w,  the  itterative  scheme  is  based  on  solving 
for  aj.  From  Figure  C-5 

3i  =  <t>0  -  aj 
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However,  j  should  be  defined  in  terms  of  the  corrected  pitch  angle  $0  cascade 
in  place  of  the  pitch  angle  <t>0;  hence,  our  Pj  is  given  by 

pj  =  <t>0  cascade  -  aj  (C.3) 

The  induced  angle  aj  can  be  expressed  in  terms  of  p j ;  Burril  (1944) 
gives  the  following  expression  for  aj 


360  Ke  T  tan(pj-p) 

aj  =  - - - sin  pj  tan  (Pj-P)  1 -  (1-KPj) 

ft  Kg5  ^scs  _  tan  pj 


(C  .4) 


The  factor  Ke  accounts  for  a  spanwise  elliptic  distribution.  Langan  and 
Wang  (1972)  show  comp ari sons  of  experimental  and  theoretical  lift  distributions; 
the  elliptic  distribution  is  good  for  r/R  <  0.95  but  for  greater  values  deviates 
noticably  from  the  experimental  results.  For  r/R  >  0.95  the  lift  distribution 
is  greater  than  or  equal  to  the  r/R  =  0.95  value,  and  the  value  at  r/R  =  0.95  is 
a  good  approximation  to  the  average  lift  distribution.  Moreover,  for  a  ducted 
propeller  there  is  even  less  likelyhood  that  the  distribution  would  decrease  to 
zero  like  an  elliptic  distribution  at  the  tip.  Since  we  will  be  using  a  nuneri- 
cal  integration  scheme  to  determine  the  thrust  and  since  the  lift  distribution 
over  the  tip  region  is  better  represented  by  the  lift  distribution  at  r/R  = 

0.95,  the  value  of  Ke  for  r/R  =  0.95  is  used  for  the  value  of  Ke  at  r/R  =  1.  Ke 
is  a  function  of  the  angle  e; 

e  =  Pj  +  (pj  -  p)  (C  .5) 


Kpj  is  the  Goldstein  correction  factor  for  a  three  bladed  propeller.  Ks  is 
a  slope  correction  factor  and  Kgs  is  an  additional  cascade  correction  factor. 
These  factors  are  to  be  found  in  Burril  (1944). 

Equations  (C.l)  thru  (C.5)  form  the  basis  for  an  itterative  scheme  to 
determine  aj,  pj,  and  <p0  cascade.  To  start  the  scheme  assume  pj  =  p  and  compute 
$0  cascade  use  (C.2)  to  get  a  starting  value  for  aj  or  use  a  value  of  aj  between 
-5’  and  15’.  In  either  case  use  (C.3)  to  compute  Pj  and  (C.4)  to  compute  aj 
out.  If  this  value  does  not  agree  with  the  starting  value  of  aj,  use  it  for  a 
new  value  of  ajjn  (i.e.  aj  going  into  the  itteration)  provided  it  lies  in  the 
range  -5°  to  15J;  otherwise,  use  another  angle  out  of  this  range.  Compare  the 
new  “lout  with  the  «iin.  A  good  common  sense  choice  of  a^n  quickly  develops, 
and  convergence  occurs  rapidly.  There  is  no  need  to  change  <f>0  cascade  or  the 
K's  for  each  new  pj;  simply  wait  for  convergence  with  the  wrong  set  of  K's  and 
than  adjust  the  K's  and  4>0  cas.  for  the  new  Pj.  After  the  first  few  times 
through  there  is  little  change  in  these  factors. 
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Table  C-5  shows  the  final  results  for  $0  cascade  and  Table  C-6  presents  the 
final  iterative  results  for  aj . 


TABLE  C-5  CORRECTFP  PITCH  ANCLE 


r/R 

0 

3j(  deg) 

Kqao 

ao  act 

4>0C  deg) 

$0cas.( deq) 

0.4 

1.01 

14.98 

0.43 

3.85 

18.89 

17.24 

0.5 

0.88 

14.30 

0.36 

3.33 

18.15 

16.96 

0.6 

0.77 

14.29 

0.28 

2.63 

17.87 

17.14 

0.7 

0.72 

14.12 

0.24 

1.55 

17.03 

16.66 

0.8 

0.65 

14.68 

0.20 

1.03 

17.22 

17.02 

0.9 

0.59 

15.79 

0.15 

0.84 

17.49 

17.36 

1.0 

0.54 

16.95 

0.13 

0.65 

18.12 

18.03 

Equation  C-2  is  used  to  compute  a  from  the  values  of  c  in  Table  C-4.  D  = 
8.50  ft.  The  hydrodynamics  pitch  anqle  Sj  comes  out  of  the  iteration;  the 
values  in  the  table  are  the  final  values.  Kqa0  is  from  Rurril  (1944),  aoact  is 
from  Table  C-2  and  <j>0  from  Table  C-4. 

The  final  values  of  the  factors  which  enter  into  the  itterative  calculation 
of  Eouation  C-4  are  tabulated  in  Table  C-6. 

The  theretical  lift  coefficient  Zvai  needs  an  experimental  correction  Ks 
and  a  cascade  correction;  the  corrected  lift  coefficient  C|_ '  is  given  by: 

Cl'  =  2u  Ks  Kgs  aj  (rad) 

These  corrections  can  be  thought  of  as  corrections  to  the  slope,  thus  the 
subscript  s. 

There  are  three  components  to  the  drag  coefficient.  Burril  (1944)  qives 
the  minimum  drag  component  in  terms  of  the  following  emperical  formula: 
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TABLF  C-fi  INDUCED  ANGLE  OF  ATTACK 


Cn  .  =  0.0056  +  0.01  t/c  +  0.10  ( t/c)2  +  «2 

Hm  n 

where  the  constant  K3  can  be  determined  from  Figure  C-9.  The  values  of  «2  in 
Table  C-7  were  obtained  through  linear  i nterpulation.  A  second  factor 

\Cnn  =  n  tan  (pj-p)  CL 
180 

is  the  induced  drag  component.  The  third  component  is  a  corrponent  depending  on 
how  close  to  optimal  lift  the  section  is  operating. 

*cD  -  K3  <C1l  .  clopt)2 

This  component  would  be  zero  at  optimal  lift  Ci  ;  qiven  together  with  K3 
Burril  as  a  function  of  maximum  camber  and  thickness°to  chord  ratio.  There  is  a 
correction  to  lift  due  to  the  drag;  Burril  gives  this  correction  as 

6  Cl  =  KqS  CD  tan  (Pi  '  P) 

~r 

Table  C-7  presents  the  values  of  all  these  lift  and  drag  components 
together  wifh  the  final  lift  and  drag  coefficients,  aj  in  radian  is  used  in 
computing  C  L;  otherwise,  the  table  is  straight  forward. 

Now  the  incremental  lift  dL  =  V2  P  V2  c  CL  i s  perpendicular  to  the  line  of 
action  of  the  resolent  velocity  V j ,  and  the  incremental  drag  dD  =  1/2  pV2  c  Cp  has 
the  same  direction.  Figure  C-5.  The  incremental 


r igure  L-4  LIFT  AND  DRAG  ORIENTATION 
thrust  dT  and  the  tangental  force  dQ/r  are  related  to  dL  and  dD  by: 
dT  =  dL  cosPi  -  dD  sinPj 

dQ  =  dL  sinPj  +  dO  cosPj 
r 
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TABLE  C-7  LIFT  AND  DRAB  COEFFICIENTS 


The  total  propeller  thrust  is 


T  =  Z 


tip 

f  dT 

root 


where  Z  is  the  number  of  blades, 
r  or  in  terms  of  the  discrete  values  in  this  report 

T  =  3  *  0.1R  l  dT 

Similarily,  the  total  torque  is  given  by 

t  i  p 

Q  =  Z  jv  dQ, 
root 


or 

Q  =  3  *  O.lR  l  dQ 

Results  based  on  these  formulae  are  presented  in  Table  C-8.  These  results 
are  computed  for  a  water  density  p  =  1.94  lb  sec  Vft^.  The  values  of  dT,  dQ/r, 
and  dQ  for  r/R  =  1.0  have  been  cut  in  half  to  reflect  that  the  interval  over 
which  these  increments  act  is  half  the  size  of  the  other  intervals.  The  inter¬ 
val  t  around  r/R  =  0.4  enters  fully  into  the  calculation;  hcwever,  the  blade 
for  r/R<0.35  to  the  hub  has  been  neglected.  Note  that  '.he  profiler  hub  has  a 
1.0625  foot  radius  which  corresponds  to  an  r/R  =  0.25.  By  neglecting  the  por¬ 
tion  of  the  blade  inboard  of  r/R  =  0.35  we  are  then  neglecting  0.425'  =  5.1"  of 
the  blade.  Within  this  segment  the  hydrodynamics  is  not  as  straight  forward. 

As  can  be  seen  from  the  results  at  r/R  =  0.4,  the  loading  is  small,  and  the 
results  will  not  change  significantly. 
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TABLE  C-8  PROPELLEP  FORCES 


The  velocity  distribution  is  obtained  usinq  the  Goldstein  (195?)  approxima¬ 
tion.  In  this  method  the  velocity  on  the  surface  is  given  by 

V  .  eCo  M  +  En  1/n  l/i  .  r,2^  sin  ( *  +  *  +  B)  -  CL  cos  ( *  +  e  +  3) 


H  =  eL0  (1  +  £ 1 )  1  l/p 
( +  sin2$>) 

+  £k 

2neCo 


A  -  £l 

\  a° 


(C.6) 


where  <(>  and  \|>  are  transform  plane  coordinates  given  by 
x  =  1/2  C  (1  +  cos  <j>) 
y  =  V?  C  U)  sin  <t> 

The  angle  e  ( <j>)  is  given  explicitly  in  this  approximation  by 


e  U)  =  (l/2n)  P  *  *  (t)  cot  1/2  ( 4>-t)  dt. 


(C.7) 


Cn  is  a  constant  for  each  section,  and  it  is  given  by 


C0  =  1 /?  (M/2)  Ip  *  (*)  d<J> 


(C.8) 


If  6  satisfies  the  condition 


1  .  Cl2V£  sin  (e  +  8)  -  CL  cos  (e  +  6)  +  C|_ 
an2/  an  2n  eCo 


(C.9) 


the  Kutta-Joukowski  condition  is  satisfied.  Ci  is  the  left  coefficient  from 
Table  07,  and  ap  is  the  left  slope  corrected  for  viscosity;  a0  =  2n  Ks  where  K 
i  s  from  Table  06. 

8  is  determined  through  an  itteratine  solution  of  Equation  (C.9).  C0,  the 

angle  e,  and  its  derivative  e'  are  computed  through  the  methods  suggested  by 
Thwaites  (1960).  The  results  for  sections  at  r/R  =  0.5,  0.7,  and  0.9  are  pre¬ 
sented  in  Table  C-9  through  Table  C-ll. 


TABLE  C-9  SECTION  BELOCITY  DISTRIBUTION  AND  PRESSURE  COEFFICIENTS  FOR  r/R 

C,  =  0.146 


TABLE  C-10  SECTION  BELOCITY  DISTRIBUTION  AND  PRESSURE  COEFFICIENTS  FOR  r/R 
Ci  =  0.166 


1.000  4.476 


From  the  Bernulli  equation 

Patm  +  V2  pV2  +  yh  =  patm  +  P]  +  V2  p  v2  +  yh 

where  gi  is  the  local  gage  pressure  on  the  foil  minus  the  static  head,  y  =  62.4 
lbs/ft  is  the  specific  weight  of  fresh  water,  and  pa^m  =  14.7  psi  =  2116.8  psf 
is  the  standard  atmospheric  pressure,  p  =  1.94  lb  sec  /ft  .  The  water  depth  h 
will  be  taken  as  the  depth  of  the  propeller  axis  below  the  design  water  line, 

7.5  feet,  minus  the  section  radius. 

h  =  7.5'  -  r 

The  local  pressure  p^  in  terms  of  the  local  velocity  v  is  given  by 
Pi  =  V2  P  (V2  -  v2); 

a  local  pressure  coefficient  Cp  can  be  defined  by 

Cp  =  PiA£  P  V2  =  1  -  (v/V)2  (C.10) 

Values  of  the  pressure  coefficient  for  r/R  =  0.5,  0.7,  and  0.9  may  be  found  in 
Tables  C-9  through  C-ll. 

For  cavity  free  operation  the  local  absolute  pressure  must  exceed  the  vapor 
pressure  e,  that  is 

Patm  +  Pi  +  Th  >  e 


or 


Cp  >  e  ~  Patm  ~  Y*1 
V2  P  Va2 


(C  .11) 


Our  analysis  of  cavitation  is  based  on  a  sunnier  time  operating  temperature  of 
70°  F;  at  this  temperature  e  =  0.36  pse  =  51.8  psf.  Table  C-12  presents  con¬ 
ditions  that  the  pressure  coefficient  must  meet  for  cavity  free  operation.  The 
velocities  V  are  from  Table  C-3. 


TABLE  C-12  CAVITATION  CRITER 


r/R 

yh(psf) 

e-yh-pgpnlpsf) 

V{ f t/ sec) 

"  2 

V2p  v  (psf) 

- " 

Cp  > 

0.5 

335.4 

MM 

3549 

-0.676 

0.7 

282.4 

-2347 

6926 

-0.339 

0.9 

229.3 

-2294 

108.61 

11442 

-0.200 

Cavitation  occurs  at  the  face  side  of  the  leading  edge  at  r/R=0.9  to 
r/R=0.5;  although  the  actual  value  of  Cp  at  r/R=0.5  is  not  below  the  critical 
value,  it  is  very  marginal.  Cavitation  also  occurs  on  the  back  side  at  r/R=0.7. 
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APPENDIX  D:  STRUCTURAL  ANALYSIS  OF  THE  BLADE  UNDER  ICE  LOADING 


As  mentioned  previously  ice  loading  is  difficult  to  predict;  as  before,  the 
criterion  proposed  by  Graft  (1977)  is  used.  The  blade  should  be  able  to  support 
a  load  equal  to  the  crushing  strength  of  the  ice  times  the  appropriate  projected 
area.  The  rotating  blade  could  be  crushing  through  an  ice  chunk  in  the  way  of 
its  rotational  motion;  in  this  case  the  projected  area  is  in  a  radial  plane  that 
is  a  plane  containing  the  axis  of  the  shaft  and  a  radial  line.  Ice  could  be 
forced  against  the  blade  in  a  forward  or  aft  direction,  in  which  case  the  pro¬ 
jected  area  of  the  blade  would  be  the  appropriate  area  to  use;  however,  in  our 
analysis  the  developed  area  is  used  for  convienance  of  calculation.  The  substi¬ 
tution  of  developed  for  projected  area  is  slightly  conservative  in  terms  of  pre¬ 
dicted  loading  and  stresses. 

Tables  D-l  through  D-3  present  the  section  properties  for  r/R=0.4,  0.6,  and 
0.8.  The  first  eight  columns  of  each  table  correspond  to  the  respective  column 
in  Table  A-l.  Ix  is  the  moment  of  inertia  around  the  cejitroidal  axis  parallel 
to  the  chord.  The  remaining  columns  are  used  to  obtain  x,  Iy,  and  Ixy,  where  Iy 
is  the  moment  of  inertia  around  the  remaining  centroidal  axis.  y 


TE 

7  =  7f  S  x  dA 
L.E. 


I 


y 


0.05  C/A  l  x  (yu-yi) 


dA 


0.05C  l  (x-x)2  (yu-yi) 


and 


Ixy  =  0.05C  l  (x-x)  (y-y)  (yu-yi ) 

Table  D-4  presents  the  corresponding  data  for  an  NACA  0018  synme_tric  foil 
section  with  a  1.5_ft  chord.  Only  yu  is  given  since  y]=-y,,.  Also  y^  for  each 
section  is  0  as  is  y;  thus,  (yi-y)  dA  is  also  zero.  Ix  is  dependent  solely  on 

Eli- 
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TABLE  O-la  SECTION  PROPERTIES  r/R=0.4 


TABLE  D-lb  SFCTION  PROPERTIES  r/R=0.4 
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TABLE  [)-2a  SECTION  PROPERTIES  r/R=0.6 
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TABLE  0-2b  SECTION  PROPERTIES  r/R=0.6 
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TABLF  D- 3b  SECTION  PROPERTIES  r/R=0.8 


0.0083ft  Iv  =  0.447  ft  =  9273  in 


TARLF  0-4  SECTION  PROPERTIES  FOR  ICF  DEFLFCTOR  STRUT 


0.0012  ft 


Actually,  the  moments  of  Inertia  around  these  axes  are  needed  but  with  the 
section  rotated  through  an  angle  equal  to  the  pitch  angle.  If  the  new  coor¬ 
dinates  are  denoted  by  x'  and  y' . 

x'  =  X  cos  <J>  -  y  sin  ij> 

y'  =  X  cos  ^  +  y  cos  <p 

rhe  desired  moments  of  inertia  Ix<  an(j  jyi  are  expressible  in  terms  of  the  pre¬ 
viously  calculated  lx,  ly,  and  Ixy. 

Ix’  =  A  y’2  dA 

_  A  (x  sin  <fi  +  y  cos  <J> )2  dA 

=  A  x2  sin2  $  +  y2  cos2  <p  +  2xy  sin  <{>  cos  $  dA 

Ix'  =  Iy  sin2  <f>  +  Ix  cos2  <j>  +  2  Ixy  sin  ^  cos  <p 

Simi  larly , 

Iy1  =  Iy  cos2  <j>  +  Ix  sin2  <f>  -  2  Ixy  sin  4>  cos  $ 

The  loading  is  also  computed  with  the  section  at  the  pitch  angle  In  this 

position  the  force  normal  to  the  x'  -  axis  at  a  given  value  of  r/R  is: 


=  (J.L  wRcj  cos  (Jii 


where  w  is  the  loading  per  unit  area  and  c;  is  the  corresponding  chord  length. 
The  force  parallel  to  the  x'  -  axis  is: 

Xj  =0.1  wKcjl  sin  <j>^ 

The  moments  that  these  forces  exert  on  each  section  are  tabulated  in  Table 
D-5  for  an  ice_crushtng  stress  w  =  300  psi.  Ad  and  A,  are  the  outboard  pro¬ 
jected  areas;  rd  and  r.  are  the  radial  locations  of  the  centriod  of  the 
respected  areas. 

The  maximum  tensile  stress  will  occur  at  the  leading  or  trailing  edge 
depending  on  whether  the  ice  is  pressing  through  the  blade  from  the  forward  or 
aft  side  of  the  screw.  If  z  denotes  the  maximum  distance  from  the  centroid  to 
the  leading  or  trailing  edge  parallel  to  the  chord, 

Zb  =  (C/2  +  x)  sin  $ 

when  computing  the  stresss  due  to  the  moment  Mj  and 
Zi  =  (C/2  +  1c)  cos  <p 

when  dealing  with  The  stress  is  given  by: 
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Tab  1  o  D-6  qives  a  tabulation  of  the  final  stresses:  <jh  hue  to  RL  and  01 
due  to  Mi. 

TABLE  D-6  ICE  CRUSHING  STRESSES 


r/R 

aj  (psi) 

°1  (psi) 

0.4 

30300 

12000 

0.6 

18600 

5600 

0.8 

8800 

2900 

The  ice  deflection  strut  is  a  beam  fixed  at  both  ends.  Roark  (1954)  qives 
the  moment  as  1/8  wl  where  the  length  is  taken  as  1  =  6'  =  72  in  and  the  loadinq 
per  foot  is  300  psi  *  1.5 1  *  12"  =  6400  Ibs/in. 

Hence, 

M  =  1/8  6400  lbs/ in  *  72  in  =  57,600  in/lbs. 

The  stress  is  given  by: 

M  Yu  max  57,600  in/lb  *  0.27  ft  *  12" 

a  =  -  =  — — - 

lx  24  i n4 

or 

a  =  7780  psi 


65 


TABLE  0-5  ICE  CRUSHING  MOMENTS  ON  BLADE 


INITIAL  DISTRIBUTION  LIST 


Command  a  lit 

U.  S.  Coast  Guard 

Washington,  D.C.  20563 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

Assistant  Librarian 
Technical  Processing  Division 
U.  S.  Naval  Academy 
Annapolis,  Maryland  21302 

Academic  Dean 
U.  S.  Naval  Academy 
Annapolis,  Maryland  21402 

Director  of  Research 
U.  A.  Naval  Academy 
Annapolis,  Maryland  21402 

Division  Director 

Division  of  Engineering  and  Weapons 
U.  S.  Naval  Academy 
Annapolis,  Maryland  21402 

Department  Chairman 

Naval  Systems  Engineering  Department 
U.  S.  Naval  Acadeny 
Annapolis,  Maryland 


Professor  T.  J.  Langan 
U.  S.  Naval  Academy 
Annapolis,  Maryland  21402 


